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The Burnett coefficient B is investigated for transport in one-dimensional quantum many-body 
systems. Extensive numerical computations in spin- 1/2 chains suggest a linear growth with time, 
B{t) ~ t, for non-integrable chains exhibiting diffusive transport. For integrable spin chains in 
the metallic regime, on the other hand, we find a cubic growth with time, B{t) ~ —D^t^, with the 
proportionality constant being simply a square of the Drude weight Dm- The results are corroborated 
with additional studies in non-interacting quantum chains and in the classical limit of large-spin 
chains. 
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Introduction. — Understanding diffusion in determin- 
istic Hamiltonian systems is one of the most ubiquitous 
problems of statistieal physics Its origin is often 

traced back to microscopic Hamiltonian chaos and its 
fundamental characterization goes beyond a single num- 
ber, the so-called diffusion coefficient D. The first higher 
order correction is systematically described by the Bur- 
nett coefficient B [lH3| which may or may not exist even 
for dynamical processes where the diffusion coefficient D 
exists 0. Even though the existence and calculation of 
Burnett coefficients for various classical statistical mod- 
els, in particular Lorcntz-type gases, has been extensively 
discussed in the literature 0, essentially nothing is 
known about Burnett coefficients in quantum systems. In 
particular recent advances in understanding the diffusion 
constant in interacting quantum spin chains @] call 
for further characterization of such transport processes 
in terms of Burnett coefficients. In this Letter we do the 
first steps by numerically simulating finite-time Burnett 
coefficients for the integrable XX Z spin-1/2 chain and 
its non-integrable modification. We find generally that 
Bit) (xt for all cases which are characterized by a finite 
diffusion constant D = Ymvt^oo D{t), not only in spin- 
1/2 chains but also in more abstract models of quantum 
transport and in the classical limit of large-spin chains. 
For integrable systems with ballistic transport behavior, 
on the other hand, we find B{t) ^ —Df-^t^ with the pro- 
portionality constant being simply a square of the Drude 
weight I?in. 

Diffusion and Burnett coefficient. — Following 0, 0, 
the (quantum and time-dependent) Burnett coefficient 
B{t) may be defined by formally expanding the decay 
rate of density-density correlation functions in Fourier 
(g— )space 



FIG. 1. (color online) Burnett coefficient B{t), divided by t^ , 
for the spin-1/2 XX Z chain at A = 0.5 for L = 8, 10, . . . , 18 
and all magnetization sectors. At short times or at A = (free 
fermions), B{t)/t^ is given by the function —t/16 (dashed 
curve). In addition, a function —(0.63)^4/16 is indicated (dot- 
ted curve), see prediction ((Sj. 



Here, D{t) is the time-dependent diffusion coefficient 
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dti/(ti), /(ii) = (J(O)J(ti)) 



(2) 



p,it)^[-q^Dit) + q*Bit) + ...]p,it). 



(1) 



given as a time integral over the two-point correlation 
function of the current operator J(t) in the Heisenberg 
picture, where (•) denotes an expectation in equilibrium, 
say, thermal state. The further occurring prefactor 
is a constant and denotes the "static susceptibility' 
The time-dependent Burnett coefficient is the difference 
between two contributions, 

Bit) = hit) - l2^2it) , (3) 
where the first term I^it) is given by 

t tl t2 

hit) = - Jdh I d<2 Jdt3 {JiO)Jiti)Jit2)Jit3)) (4) 
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FIG. 2. (color online) Absolute value of the derivated Burnett 
coefficient, |dB(t)/dt|, for the spin-1/2 XXZ chain at A = 1.5 

(a) , and for its non-integrable modification with A2 = 0.5 

(b) . Numerical results (solid curves) for L — 10, 12, . . . , 18 are 
shown in a log-log plot. In addition to the short-time behavior 
(dashed curves), the long-time behavior is extrapolated in (b) 
using the observed exponential scaling with L (dotted curve). 
Inset; The contribution l4,{t) by itself increases with L. 



as a triple-time integral over the time-ordered four-point 
current autocorrelation function. The second term 



C CI C2 

l2-2{t) = ^/dii y"dt2 Jdt3 [fih 



)f{t2 - h) 



+ f{t2)f{h-h)+.f{h)f{h-t2)] 



(5) 



is a similar time-ordered integral but over products of 
two-point correlations Conveniently, this contribu- 
tion can be rewritten as l2-2{t) = 2xD{t) dtiD{ti), 
particularly unveiling the linear increase l2-2{t) oc t in 
the case of a existent diffusion constant. But, despite 
the apparent divergence of l2-2{t) in that case, the Bur- 
nett coefficient can still remain finite, as discussed in the 
following. To this end, assume for the moment that (i) 
the two-point correlation f(t) is a delta function S{t) and 
that (ii) the four-point correlation {J{0)J{ti)J{t2)J{t3)) 
can be factorized as f{t3)f{ti — 12). Then the contri- 
butions l2-2{t) and Ii{t) are identical and, as a conse- 
quence, the Burnett coefficient vanishes exactly. Giving 
up the assumption (i) by broadening the delta function 
still allows for a finite Burnett coefficient. While the as- 
sumption (ii) appears to be crucial, it may be fulfilled for 
a non-integrable model with J, in the energy eigenbasis, 
being a random Hermitian matrix and, consequently, 
being close to proportional to an identity matrix. 



Anisotropic spin-1/2 Heisenberg chain. — We are go- 
ing to investigate the transport of magnetization in the 
spin-1/2 XXZ model as a paradigmatic example of an 
interacting manj/-particle quantum system in one dimen- 
sion. The XXZ Hamiltonian is given by 



H 



L 

E 

r=l 



(qX qX I qy qV 1 A 



(6) 



where S^'"^'^ are the components of spin-1/2 operators 
at site r, L is the number of sites arranged periodically, 
L + 1 = 1, and A is the anisotropy. In this Letter we 
consider high temperatures, where (•) = 2~-^tr(»). The 
magnetization current 



T „ ( qX qy qy qX \ 

J — 2_^l*-'r "^r+l '-'r'-'r+lJ 
r=l 



(7) 



commutes with H in the non- interacting case A = 0. In 
that case (due to (J^) L/8, (J*) = 3{L^ - i)/64, 
and X = L/A), one obtains directly D{t) — t/2 and 
B{t) = — t'^/lG, which for A > remains only an ap- 
proximation at short times, in agreement with Eq. (4) of 
Ref. 0. Remarkably, at A = a series expansion of den- 
sity autocorrelations (Bessel functions [101) leads also to 
q^D{t) and ~q^B{t) as the leading terms, hence being a 
convincing consistency check of the present approach. 

In the metallic regime, < |A| < 1, the magnetization 
current is still partially conserved such that the two-point 
correlation fit) does not decay to zero but remains at a 
finite Drude weight Dm = hmt^oo/(t), < Dm < 1/8, 
recently proven by establishing positive lower bounds in 
the thermodynamic limit [l2| . This finite Drude weight 
implies the linear increase of the diffusion coefficient at 
long times, D{t) cx 4_Dmi, just as in the case of A = 0. 
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Factoring the four-point correlation at long times 
one derives the asymptotics of the Burnett coefficient as 

B{t) ^ADlt' . (8) 

In Fig. [T] we demonstrate this result by numerically sim- 
ulating B{t) for finite length L = 10, . . . , 18 using all in- 
variant subspaces (translation, magnetization) and also 
using 4th order Runge-Kutta integration Q (step size 
5t = 0.02). While Fig. [I] clearly shows for A = 0.5 a 
stronger than quadratic increase of B(t) with time, it also 
is consistent with B(t) ^ — (0.63)^t'^/16 at long times, 
e.g., 63% of the Drude weight in the case of A = 0. Re- 
markably, the exact Drude weight at A = 0.5 for finite 
L = 18 [20] is 63% [62%] of the Drude weight at A = 0, 
while the lower bound in the thermodynamic limit is 56% 

Eventually, we discuss the regime A > 1, where Drude 
weights are expected to vanish in the thermodynamic 
limit and strong evidence of magnetization diffusion has 
been found in non-equilibrium bath scenarios on the ba- 
sis of the Lindblad equation [l^, 14]. The diffusion co- 



efficient has been shown to behave as D{t > 3.0/A) 
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0.88/ A at vanishing [HI and finite momentum @. We fo- 
cus on tfie Burnett coefficient B{t) in tliis Letter. Figure 
[5] (a) sliows numerical results summarizing an order of a 
CPU-year of simulations and plotting, for convenience, 
\dB{t)/dt\ in a log- log scale. Several comments are in or- 
der: First, after the already discussed behavior at short 
times, the Burnett coefficient changes its sign, visible as 
the zero drop in Fig. [2] (a), and indicates a correction 
towards an insulator. Second, curves for L > 12 have 
converged for times after the zero drop and show at least 
the tendency to form a plateau at i ~ 4 for L — > oo, 
then indicating a linearly increasing Burnett coefficient 
B{t) oc t. And third, even though a possible plateau is 
not pronounced yet, the contributions hit) and l2^2{t) 
increase linearly with time at < ^ 4, see the inset of Fig. [2] 
(a). Notably, the Burnett coefficient turns out to be the 
small difference between both contributions, which by 
themselves diverge with L. This divergence does not ap- 
pear in the modular quantum system, discussed later. 
Finally, we show in Fig. [2] (b) numerical results for the 
modified non-integrable XX Z model H + A2 J2r '5'r 5*^+2 
with A2 = 0.5. While the overall structure of \dB{t)/dt\ 
is similar, an emerging plateau is more clearly visible for 
the accessible lengths, hence pointing towards a linearly 
increasing Burnett coefficient again. Summarizing, a lin- 
ear asymptotic scaling B{t) oc t is clearly suggested in 
either integrable or non-integrable regimes with a finite 
diffusion constant. 

Heisenberg chains in the large-spin limit. — In addition 
we present results on the classical limit of the considered 
spin chains, where we focus on the case A = 1.5 and 
A2 = 0.5 only. In that limit the magnetization current 
is a function of classical unit (angular momentum) vec- 
tors. Their dynamics we obtain by numerically integrat- 
ing the corresponding Hamiltonian equations of motion. 
Formally, the diffusion coefficient in Eq. ([2]) and the Bur- 
nett coefficient in Eq. ([3]) remain defined the same way, 
but the equilibrium average at high temperatures is now 
performed by sampling over a set of completely random 
initial configurations, see Ref. 15| for instance. While 



the required number of initial states decreases with the 
chain length L for the evaluation of the diffusion coeffi- 
cient [l^l, the situation turns out to be different for the 
evaluation of the Burnett coefficient. As in the quan- 
tum case, the Burnett coefficient is the small difference 
between two contributions, which by themselves diverge 
with L. Thus, errors due to insufficient averaging in- 
crease with L and can only be compensated by taking 
into account more and more initial configurations. Ap- 
proximately iV « 10^ initial states are already required 
for a chain of length i = 90 taking about a CPU-year 
of computation time. In Fig. [3] we summarize these re- 
sults for three different sizes L = 18, 36, 90. Apparently, 
D{t > 10) becomes independent of time for all consid- 
ered lengths L > 18. Remarkably, the quantitative value 
D « 0.65 is close to the expectation for the diffusion con- 
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FIG. 3. (color online) Classical (a) diffusion coefficient D{t), 
(b) the contribution liit), and (c) Burnett coefficient B{t) for 
the XXZ model at A = 1.5 with A2 = 0.5. Numerical results 
(solid curves) on finite length L = 18, 36, 90 are obtained by 
numerically integrating the Hamiltonian equations of motion 
(fixed step size St = 0.05) and by averaging over A'' ~ 10* 
completely random initial states. The average over only TV = 
10^ <C 10* is also indicated (symbols). In (c) results for a 
2.5 X smaller time step are shown (dotted curves). Straight 
lines (dashed curves) serve as guides to the eye. 



stant in the quantum case [1^1 • Furthermore, B{t > 10) 
is observed to increase linearly in time, at least for the 
largest two sizes L > 36. The latter observation indicates 
the linear divergence of the Burnett coefficient in diffu- 
sive classical spin chains. This is in clear agreement with 
the finite-size results in the quantum case, see Fig. [21 in- 
dicating that the underlying mechanism "survives" the 
transition to the classical limit. 

"Modular quantum system". — Due to a delicate coun- 
terbalance of the terms I4 and /2-2, the calculation of 
Burnett coefficients for many-body systems is extremely 
demanding. Thus, in order to corroborate our prediction 
that B{t) oc t for one-dimensional lattice systems with 
finite diffusion constants, we make another numerical ex- 
periment in a single-particle diffusive quantum system 
- the so-called modular quantum system \v^ . Each 
of the L local modules features an identical spectrum, 
consisting of n equidistant levels in a band with the 
width 5e. Therefore the local Hamiltonian at the po- 
sition r is given by hr = '^^^5e/n\r,^){r,ji\ in the 
one-particle basis I?", /i). The nearest-neighbor interac- 
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FIG. 4. (color online) Diffusion coefficient D{t) and Burnett 
coefficient B{t) for the modular quantum system at coupling 
(a),(c) A = 0.0005 and (b),(d) A = 1.0. (Otlier parameters: 
n = 500, 6e = 0.5.) Numerical results (solid curves) on ar- 
bitrary length L are shown for a representative translation 
subspace (momentum k — it/ 5). Results agree well with the 
theoretical predictions (non-solid curves). 



tion between two local modules at the positions r and 
r + 1 is Vr = X rUr + H.c., 



rrir 



(9) 



with the overall coupling strength A. The r- independent 
coefficients c^.^ are a single realization of complex, ran- 
dom, and uncorrelated numbers: their real and imagi- 
nary part are both chosen corresponding to a Gaussian 
distribution with the mean and the variance 1/2. Of 
particular interest is the probability for finding the par- 
ticle somewhere within the rth local module. The as- 
sociated local current is = « A nir + H.c. with a form 
very similar to Vr, e.g., almost completely random (apart 
from the translation invariance and the necessary Hermi- 
tian property). 

The modular quantum system is one of the few quan- 
tum models which have been reliably shown to ex- 
hibit diffusion with a finite diffusion constant, reading 



D^Y^it > n/Se) = 27rA^n/(5e for weak coupling [17| and 



f^ij = {t > l/X^/2n) ~ X^'iTn/2 for strong coupling 
One might expect a finite Burnett coefficient due 
to both the presence of diffusion and the random ele- 
ments of the current. For instance, because is close 
to an identity matrix, one may be tempted to factorize 
as {J{Q)J{h)J{t2)J{h)) = f{h)f{ti-t2), then allowing 
for a finite Burnett coefficient. But the latter factor- 
ization already fails when all time arguments are equal. 
In fact, (J^) = 2/(0)^ = SA^'n^, resulting from an ad- 
ditional coherent sum over the off-diagonal elements of 
J^. Instead, fulfilling the static property, we may choose 
the unbiased factorization of ( J(0) J(ti) J(t2) J(t3)) into 

2/i[f{h)f{h - h) + f[h)f{h - h) + f{h)f{h - h)], 



yielding the relation Ii{t) = 2/3/2-2 (i) between the two 
contributions to B{t). Therefore, noting that x = 1j this 
choice leads to the linearly increasing Burnett coefficient 
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B{t > t) 



(10) 



For verification, we present in Fig. 2] numerical results on 
D{t) and B{t), as obtained from the application of full 
exact diagonalization and the additional use of 4th order 
Runge-Kutta integration (step size 5t <^ t). Because the 
linear growth of the Hilbert space with L is compensated 
by translation invariance, the dimension of a momentum 
fc-subspace is only set by the level number n, chosen as 
n = 500 to ensure a sufficient number of states. Since we 
do not find a significant dependence on fc. Fig. |4] depicts 
numerical results for a single fc-subspace. The quantita- 
tive agreement with the theoretical predictions on D{t) 
and B(t) clearly demonstrates a linearly increasing Bur- 
nett coefficient despite the existence of a diffusion con- 
stant, which is in agreement with the previous results 
on spin chains, but much clearer due to computational 
simplicity of the model. 

Conclusion. — In this Letter we presented extensive nu- 
merical and theoretical investigations of Burnett coeffi- 
cients in quantum chains. We conjectured and supported 
the general statement that in the thermodynamic limit 
Burnett coefficients diverge linearly, B{t) oc t, in diffusive 
regimes with finite diffusion constants. Due to Eq. ([T|), 
this linear divergence is still consistent with diffusion but 
causes Burnett coefficients to be a relevant correction at 
arbitrary small momentum. In ballistic regimes with pos- 
itive Drude weights, on the other hand, we demonstrated 
the cubic divergence B(t) cx t^. This behavior is remark- 
ably different than for Lorentz billiard type classical sys- 
tems and and calls for a deeper theoretical analysis. 

Financial support by the grant P 1-0044 of Slovenian 
research agency is gratefully acknowledged. 
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